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ABSTRACT 

A new thermodynamically consistent damage model is proposed for the simulation of 
high-cycle fatigue crack growth. The basis for the formulation is an interfacial 
degradation law that links Fracture Mechanics and Damage Mechanics to relate the 
evolution of the damage variable, d, with the crack growth rate da/dN. The damage 
state is a function of the loading conditions (R and AG) as well as the experimentally- 
determined crack growth rates for the material. The formulation ensures that the 
experimental results can be reproduced by the analysis without the need of additional 
adjustment parameters. 

1. INTRODUCTION 

Mechanical fatigue, especially high-cycle fatigue, is a common cause of failure of 
aerospace structures. In laminated composite materials, the fatigue process involves 
several damage mechanisms that result in the degradation of the structure. One of the 
most important fatigue damage mechanisms is interlaminar damage (delamination). 
There are two basic approaches for the analysis of delamination under fatigue: Fracture 
Mechanics, which relates the fatigue crack growth rate as a function of the energy 
release rate and mode ratio; and Damage Mechanics, in which the concept of a 
cohesive zone [1,2] is used to establish damage evolution as a function of the number 
of cycles. 

In a degradation process involving high cycle-fatigue, damage evolution can be 
obtained as the sum of the damage caused by static or quasi-static loads and the damage 
that results from the cyclic loads. The damage evolution produced by cyclic loads is 
usually formulated as a function of the number of cycles and strains (or displacement 
jumps) [3-5], where a damage evolution law expressed in terms of the number of cycles 
is established a priori. However, the damage evolution law must be expressed as a 
function of several parameters that have to be adjusted through a trial-and-error 
calibration of the whole numerical model. In this paper, an alternative approach is 
proposed whereby the evolution of damage is based on linking Fracture Mechanics and 
Damage Mechanics, and relating the evolution of the damage variable, d, with the 
crack growth rate. 
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The present model is implemented by means of a user-written element in ABAQUS [6] 
by adding the damage evolution law formulated in the decohesion element previously 
developed by the authors [7,8], 

2. CONSITUTIVE MODEL FOR QUASI-STATIC LOADS 

The constitutive model used in this paper for quasi-static loads uses a bilinear 
constitutive equation that relates surface tractions, z, to displacement jumps. A, as 
shown in Figure 1. While the interface is undamaged, a high interface stiffness, K, 
ensures a stiff connection between the two adjacent layers. Interface damage initiates 
when the displacement jump norm reaches the onset displacement. A 0 , and the interface 
is considered fractured when the displacement jump norm exceeds the final 
displacement, A. The energy dissipated during the damage evolution is called E , and 
the area under the traction-displacement jump law is equal to the critical energy release 
rate G c . The critical displacement jumps A 0 and A are functions of the mode mixity and 
details of the formulation can be found in [8], 



Figure 1. Bilinear constitutive law used for quasi-static loads. 


3. CONSITUTIVE MODEL FOR HIGH CYCLE FATIGUE 

The damage evolution in a degradation process involving high-cycle fatigue can be 
considered as the sum of the damage sustained from quasi-static loads and the damage 
sustained from cyclic loads: 


^ ^static ^cyclic ^ | j 

In the framework of the Damage Mechanics, the damage evolution that results from 
cyclic loads can be formulated as a function of the number or cycles, usually as a 
function of strains (or displacement jumps) [3-5], where a damage evolution law is 
established a priori as a function of the number of cycles. The damage evolution law is 
a function of several parameters that have to be adjusted calibrating the numerical 
model with experimental results, usually by trial and error. However, in this paper, the 
evolution of the damage evolution law is formulated using Fracture Mechanics and 
creating a link between Fracture Mechanics and Damage Mechanics to relate the 
evolution of the damage variable, d , with the crack growth rate, da/dN. The evolution 
of the damage variable is related to the evolution of the crack surface as follows: 
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( 2 ) 


5d _ 5d dA d 
dN ~ dA d dN 

where A d is the damaged area, and is the growth rate of the damaged area. While 

the second term in the right hand side of equation (2) must be characterized 

dd 

experimentally, the first term can be obtained from either a Damage Mechanics 

approach or a Fracture Mechanics approach, as described in the following sections. 


3.1- Determination of using a Damage Mechanics approach 

54 , 

In the framework of the Damage Mechanics, the damage variable, d , can be defined as 
the ratio between the damage area and the total area, i.e., 


d 



(3) 


The damaged area is zero when no damage is present in the process, A d = 0 , while 
A d = A when the surface is completely damaged. Since the total area is constant, the 
variation of the damaged variable with respect to the damaged area can be written as: 


5d _ 1 

^ = A 


(4) 


3.2- Determination of using a Fracture Mechanics approach 

54 , 

In the framework of the Fracture Mechanics, the fraction of the damaged area, A d , with 
respect to the total area, A , can be written as a function of the dissipated energy: 


A 

A 



(5) 


where 5 is the fraction of the energy per surface unit dissipated during the damage 
process, i.e., the area under the cohesive law for the current damage threshold, and G c 
is the critical energy release rate (see Figure 1). 


Assuming no change between modes, G c is constant, while 5 is a function of the 

cohesive law used and the current damage threshold. Using this approach, the variation 
of the damage variable with the damaged area can be written as: 

dd_ = G^dd 

54, Ad E 
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The evolution of the energy dissipation with the damage evolution is obtained from the 
equations of the constitutive law used. Using a bilinear constitutive law: 


dS_ _ G 4 ( 4 - 4 ) 

3d c (z/(/-d) + dzf) 2 


( 8 ) 


Using equation (8), equation (7) can be written as: 

0d _ 1 _ U(7-d) + dzf) 2 
04 - A 4(4 - 4 ) 


3.3- Determination of the damaged area growth rate as a function of the number of 
cycles 


In this section, we establish the dependence of the damaged area growth rate. 


04 

dN 


on 


the surface crack growth rate, 


dA 

dN 


For a specimen with just one crack front, the crack 


growth rate is equal to the sum of the damaged surface growth rates of all elements in 
the cohesive zone. In the other regions of the specimen, there is no possibility of new 
surface generation. 


dA 

dN 


y 34 

,k » 


( 10 ) 


Using the simplification that 
equation can be written as: 



is constant over the cohesive zone, the previous 


dA 

dN 


y 34 

par 


^ -cz 34 

A dN 


( 11 ) 


where the ratio — 

A 


represents the number of areas in which the cohesive zone have 


been divided. In a finite element environment, this ratio represents the number of 
elements that span the cohesive zone. Rearranging terms in equation (11), the surface 
damage growth rate can be written as: 


04 _ A dA 
dN ~ A cz dN 


( 12 ) 


3.4- Evolution of the damage variable under cyclic loading 

Replacing equation (12) and equation (4) (based on Damage Mechanics) in equation 
(2) the evolution of the damage variable as a function of the number of cycles can be 
written as: 
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dd 1 dA 


( 13 ) 


dN A cz dN 


Likewise, using the approach based on Fracture Mechanics, equation (7), the evolution 
of the damage variable as a function of the number of cycles can be written as: 

dd_ = G^dd_dA 
dN A cz dZ dN 


where G c depends on the material used and loading mode, and — depends on the 

dZ 

cohesive law used in the formulation of the surface traction-displacement jump 
relation. 

The area of the cohesive zone can be computed intensively using Rice’s model [9]: 


A 


cz 


9 n E 3 G 
b 

32 (r’f 


(15) 


where b is the width of the specimen, G is the energy release rate, E 3 is the Young‘s 
modulus of the bulk material in the direction perpendicular to the crack plane, and z° 
is the interfacial strength. 


3.4- Crack growth rate 

dA 

The surface crack growth rate under fatigue loading, , is a load and material- 

" ~ dN 

dependent characteristic that has been widely studied, and models such as the Paris law 
have been developed to represent it. For instance, 


dA 

dN 


= C 


/ \ m 

f AG^ 

G ■ j 


(16) 


where C and m are parameters that depend on the mode ratio and must be determined 
experimentally. AG, is the cyclic variation in the energy release rate, which can be 
computed using the constitutive law of the interface 

G = \MA)dA (17) 

for the maximum and minimum displacement jump over the cycle (see Figure 2): 

AG = G -G • (18) 

max mm V ' 

Defining the reversibility factor, R , the relation between the minimum and the 
maximum displacement is: 


A min = RAT 


ax 


(19) 
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and the variation of the energy release rate can be written as: 


AG = 


KAA“ ax \Af -zf‘ 


1 + R 


{l-R) 


( 20 ) 



Figure 2. Variation of the energy release rate. 


4. RESULTS AND DISCUSION 

The present model is implemented by adding the fatigue component of damage to the 
damage evolution law of a user-written decohesion element previously developed by 
the authors [7,8]. To verify the constitutive model, several single-element tests were 
performed. The evolution of the interface traction in the constitutive equation for a high 
cycle fatigue test subjected to displacement control is shown in Figure 3. 



Figure 3. Evolution of the interface traction in the constitutive equation for a displacement 
jump controlled high cycle fatigue test. 


It can be observed that fatigue loading results in a reduction of the stiffness, and in a 
reduction of the interfacial traction and interfacial strength. The evolution of the 
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interface traction and the interfacial strength with the number of cycles is shown in 
Figure 4. The shape of the obtained curves is similar to the widely-used S-N curves 
used in the design for fatigue strength. 

Several tests were conducted to simulate the surface crack growth velocity under mode 
I loading for different levels of the energy release rate. A description of the 
experimental procedure that was simulated is reported by Asp et al. [10]. The boundary 
conditions used in the simulation are the same to those presented by Robinson et al. [4], 
The results obtained from the simulations are shown in Figure 5 against the 
experimental data. A better agreement between simulated and experimental data is 
observed when the Fracture Mechanics approach is used instead of the Damage 
Mechanics approach. 



Figure 4. Evolution of the interface traction and the maximum interface strength as a function 
of the number of cycles for a displacement jump controlled high cycle fatigue test. 



Figure 5. Comparison of the experimental data with the crack growth rate obtained from the 
numerical simulation for a Mode I DCB test. 
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5. CONCLUSIONS 


A thermodynamically consistent damage model for high-cycle fatigue delamination 
was developed. The evolution of the damage variable was derived by linking Fracture 
Mechanics and Damage Mechanics to relate damage evolution to crack growth rates. 
The damage evolution laws for cyclic fatigue were combined with the law of damage 
evolution for quasi-static loads within a decohesion element previously developed by 
the authors. The model was validated with single-element numerical tests as well as a 
simulation of a DCB problem. The model was able to reproduce the test data without 
the need of additional adjustment parameters that are typically used in other fatigue 
growth models. 
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